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Abstract. We derive semiclassical trace formulae including Gutzwiller's trace formula using 
coherent states. This formulation has several advantages over the usual coordinate-space for- 
mulation. Using a coherent-state basis makes it immediately obvious that classical periodic 
orbits make separate contributions to the trace of the quantum-mechanical time evolution 
operator. In addition, our approach is manifestly canonically invariant at all stages, and 
leads to the simplest possible derivation of Gutzwiller's formula. 
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1 Introduction 

Gutzwiller's trace formula jj], ||, || relates the density of states of a quantum system 
to periodic orbits of the corresponding classical Hamiltonian. In its usual form it 
applies to systems with a discrete quantum spectrum, and isolated, unstable periodic 
orbits. In the general case it is not an exact relation. The density of states p(E) may 
be expressed in terms of the trace of the quantum- mechanical time evolution operator 

— exp(iEt/h)Tr[U t ] (1) 

where E n are the eigenvalues and Ut = exp(— ifft/Ti) is the quantum-mechanical 
time evolution operator. The idea is to approximate the trace by a formula valid in 
the semiclassical {% — ► 0) limit. Gutzwiller has also provided a discussion of some 
precursors of his expression in the mathematical literature ||, including an exact 
formula due to Selberg (discussed in which applies to the Schrodinger equation 
on closed manifolds of constant negative curvature. Gutzwiller's formula has been 
used as a basis for many investigations in semiclassical quantum mechanics, with 
applications in atomic j^], molecular and mesoscopic physics Hj. We will discuss 
some of its extensions in due course. It is certainly a powerful and fundamental 
element in semiclassical quantum theory, and it is desirable to have a clear and direct 
derivation. 

The original derivation by Gutzwiller is undeniably far more complicated than 
the final result, and it also has the disadvantage that it uses the coordinate rep- 
resentation of the quantum-mechanical time evolution operator, whereas the final 
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answer is expressed in a canonically invariant form. The approach described here 
uses a coherent-state basis to evaluate the trace. The coherent states \a) will be 
defined later. They are labeled by a phase-space point a = (q,p) — q and p are 
d-dimensional coordinate and momentum vectors respectively — and may be thought 
of as wave packets positioned at the phase-space point a. The trace of any operator 
may be expressed as an integral over coherent states. For example, the trace of U t 
may be written as 

^]=/(Sf^w- (2) 

The semiclassical approximation for the state resulting from the action of the time 
evolution operator Ut on the coherent state \a) is clearly expected to be a wave packet 
with phase-space coordinates a t , which are the result of mapping the phase-point a 
for time t under the action of Hamilton's equations. The contribution to the integral 
(0) from the phase point a is clearly negligible, unless a is very close to a periodic 
orbit of period close to t. The fact that the density of states may be expressed in 
terms of periodic orbits is thus immediately obvious when the coherent state basis is 
used. This point was made some time ago in a paper by one of us 0. 

It is less obvious how the precise form of Gutzwiller's trace formula may be ob- 
tained using coherent states. A derivation based upon coherent states has previously 
been given by Combescure et al [0 , but we will argue shortly that our method has the 
advantage of being canonically invariant, as well as being considerably simpler. Our 
calculation builds upon the semiclassical evolution of coherent states as described in 
[pi . Translation of wave packets in phase space is effected by so-called Weyl-Heisenberg 
operators, T(a). Wave packets are not only moved around in phase space but, in 
general, also deformed under quantum evolution. This deformation is described by a 
symplectic transformation represented by a matrix S. A representation of symplectic 
transformations on the quantum-mechanical Hilbert space is achieved by metaplectic 
operators iZ(S) (they actually form a double covering of the symplectic group) . In Jo) 
the metaplectic operators are expressed in terms of coordinate-space matrix elements, 
and a similar coordinate space representation was used by Combescure et al. Our cal- 
culation uses the Weyl representation of metaplectic operators, partially developed 
in ||lol| , in which an operator A is represented as an integral over Weyl-Heisenberg 
operators T(a), with weight A(a). This representation gives a canonically invariant 
formulation, without reference to any coordinate space representations. 

Our derivation of Gutzwiller's formula is based upon a simple formula for the trace 
of a metaplectic operator 

Tv[ n { s)] = *0^B (3) 

|det(S-/)| 

(where v is an integer index). It uses the fact that the periodic-orbit contributions 
to Tr[C/t] are proportional to exp(iR p /K)Tr[R-(M p )] where R p is Hamilton's action 
along the periodic orbit and M p is the corresponding stability matrix. The use of this 
expression leads to a derivation which is arguably the simplest possible, since none of 
the intermediate stages involve expressions which are significantly more complicated 
than the final result. All of the stages are canonically invariant. 

An advantage of finding a simplified derivation is that extensions and variations 
of the original expression are easily obtained. Gutzwiller's trace formula has been 
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extended to give an expression for the density of states of integrable systems 
It has also been extended by one of us to yield information about diagonal |l2| and 
off-diagonal |?J matrix elements of operators. We will also give very direct derivations 
of the latter two formulae, and discuss the integrable case briefly in the concluding 
section. 

The remainder of this paper is organized as follows. Section |2| deals with some 
definitions and notation. Section |^ discusses a semiclassical approximation using 
the Weyl representation of metaplectic operators. This representation of metaplectic 
operators has not previously been developed in detail. It allows us to work entirely 
in phase space, thus simplifying the discussion considerably. In section ^ we show 
how to derive a trace formula for the density of eigenphases of a Floquet matrix 
corresponding to a chaotic quantum map. While this example has limited physical 



applications, it serves as a simple illustrative example. In section 4.2 we derive the 
Gutzwiller trace formula for chaotic flows. The discussion is parallel to that of section 
[|, with the complication that the marginal direction along the flow has to be treated 
separately. Section [4.3| considers more general trace formulae for chaotic systems. 
Section || contains a discussion of the results. 

2 Definitions and notation 

2. 1 Traces 

We consider a closed quantum system described by a time-independent Hamiltonian 
H, for which the solution of the stationary Schrodinger equation gives rise to discrete 
energy levels Ej and eigenfunctions H\ipj) — Ej \ipj). The density of states is 
defined as 

p(E) = J25(E-Ej). (4) 

3 

In the following it will be shown how to obtain semiclassical estimates of this and 
more general densities such as 

PA (E) = J2(^\ A \^) HE - E^ , (5) 

3 

S A {E,huj) = ^\{i>j\A\iPi)\ 2 5(hu-E k +Ej)5(E-Ej). (6) 

j'.J 

Here, A is an observable with a classical limit. The three densities p(E), Pa(E) and 
Sa{E,Tiu>) may be written as traces over the quantum-mechanical time evolution 
operator 

U t ee exp( - iHt/ti) . (7) 



One has 



/OO J i 
— exp(iEt/h) Tr [U t ] , (8) 

f°° dt 

Pa(E) = J —exp(iEt/n)Tr[AU t ], (9) 

f°° dt f°° ds .... 

Sa(E, huj) = / — exp(iEt/h) — exp(hjs)Tr[AU a AUt- a ]. (10) 
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In the remainder of this article, semiclassical expressions for the traces occurring in 
Eqs. (|) to @ will be sought. 

We will also consider explicitly time-dependent systems, subject to a driving force 
with period T. In this case the density of eigenphases 9j of the Floquet operator Ut 
will be calculated. The Floquet operator maps solutions \ipt) of the time-dependent 
Schrodinger equation according to \ipt+kT) — Uj, \ipt) where k is an integer. Ut is a 
unitary operator and the density of its eigenphases is defined as 

oo N 
I — — 00 j — 1 

This may be written as 

oo N 

p(*) = i E E^ n{6 ~ 6j) ( 12 ) 

n— — oo j — 1 

using the Poisson summation formula. Thus 

1 °° AT -. 

m = 2^ E Tr Pt\ = ^ + ^ E Tr fa] j* . (13) 

71 — — 00 n>0 

In this case, we seek a semiclassical expression for Tr [U$] for n > 0. Eq. (||) may 
be viewed as a discrete analogue of Eq. (||). 



2.2 Classical dynamics 



We consider a possibly time-dependent Hamiltonian flow defined by the Hamilton 
function H(q,p; t) with d degrees of freedom. Here a = (q,p) = (qi, ■ ■ ■ ,qd,Pi, ■ ■ ■ ,Pd) 
are phase-space coordinates, a t is a family of phase-space points representing a tra- 
jectory labeled by time t. The flow is governed by Hamilton's equations 



~ r dH 
ot t = Jt. — 
da t 



with 



J = 



I 

—I 6 



(14) 



(15) 



The separation 8ot t of two nearby trajectories is given by the stability matrix M t 

Sat = Mt Sao ■ (16) 
The stability matrix M t obeys the differential equation 

JK t M t , (17) 

where we have defined the matrix K t with matrix elements 

d 2 H 



dt 



{Kt}, 



dctidctj 



a = at 



(18) 
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3 Semiclassical dynamics in phase space 

We will make use of a representation of quantum dynamics in phase space and will 
evaluate traces such as the trace Tr [U t ] of the time evolution operator U t in a coherent- 
state basis. 



3.1 Coherent states and Weyl-Heisenberg operators 

Following the notation established in H we denote coherent states by 

\a Q ) = f(ao) |0) (19) 

where the state |0) is the ground state of a harmonic oscillator and the Weyl-Heisenberg 
operator T(c*o) is given by 

t (a„) = e-K^o A «)/» (20) 

where a = (q,p). Throughout we define olq A a.\ = <7 -Pi — 9i-Po; where <ZP denotes 
the usual scalar product. Weyl-Heisenberg operators have the following properties 
(all of which may be verified using the Baker-Cambbell-Hausdorff formula |l3| ) . 

1. Translation 

= a + a . (21) 

2. The Weyl-Heisenberg operators satisfy the following concatenation property 

f (a )f (ai) - e - ^ A «i]/ 2?i f (a + ai) . (22) 

From this it follows that T _1 (ao) = ao) = T^(ao). 

3. Equation of motion 

ih^f(a t ) = f (a*) (|[a t A a 4 ] + [a t A a]) . (23) 

The Weyl-Heisenberg operators translate wave packets in phase space, by virtue of 
(pi]). It is expected that wave packets are not only moved around but also deformed 
under quantum evolution. In order to represent deformations of wave packets one 
needs a representation of symplectic transformations S on the Hilbert space of the 
system. 

3.2 Representation of symplectic transformations 

We seek operators which provide a representation of symplectic matrices S on the 
Hilbert space of the system. These are termed metaplectic operators, and are defined 
by the requirement that the operators 7Z(S) which which represent the symplectic 
matrices S are unitary matrices, satisfying 

7^( 1 S)f(a)7?.~ 1 (5) =f(,5*a). (24) 

Usually |^|, |^] this representation is derived in configuration space and explicitly de- 
fined by prescribing the configuration-space matrix elements (q'\7l(S)\q) oiTZ(S). In 
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the mathematical literature so-called Fock-space representations of symplectic trans- 
formations are often used. These go back to and are defined in [^4). Such repre- 
sentations introduce additional undesirable complications in the present context. In 
the following we use the Weyl representation jl5| of the metaplectic operators TZ(S), 
in which an operator is expressed in terms of an integral over the Weyl-Heisenberg 
operators. This representation is simpler, and also more natural because it uses phase 
space throughout, avoiding making an arbitrary choice of coordinates. The only place 
where we are aware of the Weyl representation of the metaplectic operators having 
been constructed previously is in jnj (the results differ slightly: the requirement 
of unitarity was not imposed, the expansion was in terms of magnetic translations 
rather than Weyl-Heisenberg translations, and a phase factor was included which 
always evaluates to unity). 
The operator 

is the metaplectic representative of the d-dimensional symplectic matrix S when the 
symmetric matrix A is given by 



A= ±J(S + I){S - I)- 1 . (26) 

The integer index v will be discussed further in due course, but we remark that 
when all of the eigenvalues of A are large and positive, setting v = d makes H(S) 
approximate the identity operator, /. Following the approach of reference JTo[ , using 
( p2[ ) it is easily verified that this expression satisfies (|24|), and that the following 
properties hold: 

1. The operator TZ(S) represents symplectic transformations on the quantum- 
mechanical Hilbert space, in the sense that it satisfies (0), and also 

&(5)a£ -1 (5) = fiMa. (27) 

2. The operator 1Z(S) is unitary 

(28) 

3. The operator 7t(S) obeys the following law of concatenation 

n(s 1 )n(§ 2 ) = ±n{s 1 s 2 ). (29) 

The origin of the sign will be discussed in due course. 



Thus ( p5[ ) constitutes the desired quantum-mechanical representation of symplectic 
transformations. Furthermore, the following properties will be used below: 

4. The representation of a symplectic matrix with an infinitesimal generator, S = 
I—eHJ (with H symmetric) is given by 

K(S) = I - ^a.JHJa + 0{e 2 ) , (30) 



compare |Tj| [Ti| . 
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5. Using the infinitesimal representation (^0|), it was shown in || that the operator 
Tt(St) obeys the equation of motion 

A A Q 

m-7t(5 t ) = -^a.J^-S^aliiSt) . (31) 

6. Finally, the trace of 7^.(5) is 

Tr[7^)]= ; XP(W2) . (32) 
V|det(5-7)| 

This is immediately obvious from (p5[) by noting that Tr [7"(y)] = (2nh) d S(y). 



Obviously, the representation (|25|) breaks down whenever S has eigenvalues unity. 
It will be seen below that in a semiclassical description of wave-packet dynamics, S is 
just the stability matrix of the classical nearby-orbit problem. In the case of chaotic 
Hamiltonian flows, this stability matrix always has at least two eigenvalues unity, 
corresponding to the direction along the flow (parallel direction). In this case it is 
possible, as will be shown below, to separate the metaplectic representations corre- 
sponding to motion in the parallel and perpendicular directions. For the metaplectics 
in the parallel direction, the following representation will be useful: the singular 2x2 
symplectic matrix 

' 1 b 



S=[ o i ) (33) 



is represented by the metaplectic operator 



1 C°° 

ti(S) = -== / d Xl e lx *l 2hb f{x) (34) 



with x — (xi,0) (c.f. WA 10|). An infinitesimal negative imaginary part added 
to the real number b makes this integral convergent. Note that the phase of the 
prefactor increases by tt/2 when b goes negative. This representation of a subgroup 
of symplectic transformations satisfies the above properties 1-6. 



3. 3 Topology of the symplectic group 

We must now consider the topology of the group of symplectic matrices and of the 
corresponding metaplectic operators. In particular, we must explain the fact that the 
metaplectic operators form a 'double covering' of the symplectic matrices, since this 
property gives contributions to the Maslov indices JlJ] which occur in the Gutzwiller 
trace formula. 

We commence by discussing some results which are already known || . We consider 
only the case of 2 x 2 matrices; the results generalise to higher dimensions. A general 
symplectic matrix may be written in the form 

S = S T S R (35) 

where St and Sr are, respectively, a symmetric matrix satisfying det(Sr) = 1, and 
an orthogonal matrix. In the case of 2 x 2 matrices, the matrix St belongs to a 
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two-parameter family with the topology of the plane, and Sr is a rotation matrix 
parametrised by an angle 9, so that the space of the orthogonal matrices therefore has 
a circular topology. The space of 2 x 2 symplectic matrices therefore has the topology 
of the Cartesian product of a plane and a circle. The symplectic matrices can also 
be labeled by points in a space of three real variables, as a periodic function of the 
parameter 9. It will be helpful to bear both representations in mind. We will term 
these two regions of parameter space the 'fundamental domain' and the 'extended 
domain', respectively. The extended domain consists of a set of 'cells', which are 
replicas of the fundamental domain generated by adding an integer multiple of 2ir to 
9. 

The metaplectic operators TZ(S) may be represented as analytic functions of the 
underlying symplectic matrices S, but this analytic representation could be multi- 
valued. We discuss the simplest case, in which the matrix St is the identity matrix, 
so that S is a 2 x 2 rotation matrix, Sr(6). The metaplectic operator representing 
a rotation in phase space is clearly generated by the action of a harmonic oscillator 
H(x,p) = \{x 2 +p 2 )- The corresponding metaplectic operator may be written either 
in the form (p5|), or alternatively 

n(S K (6)) = exp(^(x 2 +p 2 )9^j . (36) 

By expanding states in a basis of harmonic oscillator eigenfunctions, we see that in 
the extended domain, this operator clearly satisfies the concatenation rule in the form 
R(Sr(0i+02)) =n(S K (0 1 ))n(S K (d 2 )),bnt'R(S K (2Tr)) =-I. We therefore see that 
in the fundamental domain, the metaplectic representing rotations is a double cover- 
ing, in the sense that a product of two metaplectics is equal to another metaplectic 
only up to a sign. Specifically, for 9 in the fundamental interval [0, 2tt), we have 



K(S r (0i))1I(Sr{02)) = (-l) s K(S R ((9i + 9 2 )mod27r)) , s = intf J ^ ± ) • (37) 

This property is general, in that the metaplectic operators do form a double covering 
of the symplectic operators in the fundamental domain, although they can be repre- 
sented as a single valued analytic function in the extended domain. This assertion can 
be checked by means of the following argument. Firstly, using the properties of the 
Weyl-Heisenberg operators, it is verified that ( ^H|) satisfies (f29|), with a factor (— l) s 
appearing rather than some more general phase factor. Continuity considerations 
show that if ([57|) holds when S = Sr, it must also hold in the more general case, 
implying a double valued rather than single valued function throughout the funda- 
mental domain. The fundamental domain could be delimited by the planes 9 = and 
9 = 2ir in the representation ([55]), although there is an arbitrariness in the choice of 
the boundaries, and we will discuss an alternative choice shortly. 

We now discuss the index v appearing in our general representation of the meta- 
plectic operators, defined by (|25|). This is an explicit analytic representation, which is 
defined everywhere except upon the manifold of codimension 1 where det(S — J) = 0. 
(Alternative representations, such as ([54]), are available on this manifold.) The index 
v is chosen so that the resulting representation of the metaplectic operators is con- 
tinuous across manifolds where det(S — I) = 0. It changes by ±1 on each crossing. 
The manifolds where det(S — I) = are of two types; the function det(S — I) is 
either increasing or decreasing as 9 increases. These form an alternating sequence as 
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9 increases (although they may coalesce to a double zero, as happens when St = !)■ 
For our purposes it will be convenient to define the domain boundaries as being one 
of these two sets of surfaces satisfying the equation det(5 — J) = 0. There are two 
possibilities for the change of the index upon crossing to a successive cell. Either the 
changes of index have the same sign on each of the two manifolds where det(S-I) = 0, 
or they have different signs. The fact that the metaplectic operators change sign when 
6 increases by 2n implies that the former possibility is realised. The index v therefore 
changes by ±2 when crossing to an equivalent point in an adjacent cell. 



3-4 Semiclassical dynamics 

In Ref. H the following ansatz for the quantum-mechanical time evolution operator 
is discussed for the evolution of the coherent state la) 



U t \ ao )~e^/ n f(cc t )ii(S t )\0). (38) 



It is not possible to write a corresponding semiclassical ansatz for the operator U t 
itself, in a form which is independent of the label of the coherent state upon which it 
acts. Nevertheless, ( |38| ) is sufficient for our purposes. The phase 74 and the matrix 
St are to be determined by substitution into the Schrodinger equation 

ih^-\a t ) =H\a t ). (39) 
at 

The lhs of © gives 

-W 7t/ *t(a t )ft(S t )|0) 

+jrt/ n f(a f ) (i [a t A a t ] + [a* A a}) H(S t ) |0) 

^le^/ h f( at )n(S t )aJ^Sr 1 an(S t )\0). (40) 
z at 

The rhs of ( J39| ) is evaluated using 

H\a t ) = jYt/ n (41) 
x {H(aut)f(a t )n(S t )\0)+f(a t )—aTZ(S t )\0) 

+ \f{a t )a..J K t aR{S t )\Q) + • 

with K t given by Eq. (|l8|) . This expression can be obtained from the Weyl quantisa- 
tion scheme. One thus obtains the following set of equations, 

jt = i[QiAd t ]-ff(a ( ), (42) 
~ dH 

at = J^~, 43 

OOLt 

±St = JKtSt- (44) 
From ( p2| ) it follows that 74 is related to Hamilton's action R t 

It = Rt - + kPo-Qo > Rt = Jp-dq-Hdt (45) 
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Eqs. ( p3| ) are just Hamilton's equations, (14). By comparison of (0) and Q it 
follows that St is the stability matrix Mt- The final result is 

|ot) - e 1 ^* ~ Pt-<Zt/ 2 + Po-Qo/ 2 )/* t( at ) K(M t ) |0) . (46) 

The propagation of a state |ao) is thus envisaged as follows: \o.q) is transported to 
the origin of phase space by T^cto), there it is deformed by TZ(M t ) and finally the 
state is transported to \a t ) by T(a t ). 

The evolution of the wave packet should be continuous. The metaplectic operators 
can be presented as an analytic function of the symplectic matrices in the extended 
domain. If wc follow the evolution of the symplectic matrix as the trajectory evolves, 
the metaplectic operator evolves continuously, and expression (^) remains valid. The 
index v changes by ±1 each time M t crosses the a manifold where det(M t — I) = 0. 
These changes of index are the origin of the Maslov indices Eq] which appear in the 
Gutzwiller trace formula. 

4 Trace formulae for classically chaotic (hyperbolic) systems 

4-1 Gutzwiller 's trace formula for chaotic maps 

In this section, we consider a time-dependent Hamiltonian flow with d — 1 degree 
of freedom, defined by the Hamilton function H(q,p;t), such as a driven pendulum 
with a periodic driving force of period T. We define a stroboscopic section at t = nT, 
with n = 0,1,2,.... The dynamics on the stroboscopic section, generated by the 
flow H(q,p;t), is given by a Poincare map acting on the stroboscopic section. The 
quantum evolution operator of this Hamiltonian, U t evaluated at the stroboscopic 
time t — T is termed the Floquet operator. We will obtain an expression for the trace 
of powers of the Floquet operator U n T — [Ut] 11 , and use this to give a trace formula 
for the density of eigenphases. We may write 

Tr[L> t ]= / ^(a\& t \a). (47) 

Now 

(a\U t \a) ~ exp[i(R t - Pt q t /2 + Po qo/2)/h} (0\t*(a) f(a t ) K{M t )\0) . (48) 

(Here we consider TZ(M t ) in the extended domain.) We assume that the classical 
dynamics is hyperbolic. In this case, within the semiclassical approximation, the 
essential contributions to this trace will come from isolated, unstable periodic points 
ct p and we write (compare Fig. |l|) 

a = ct p + 5a , a t = a p + M p 8a , (49) 

where t = nT, and M p is the stability matrix evaluated at the periodic point. Using 
~22|) and ( |24| ) one has for the matrix element in the vicinity of a p 

(OIT+Op + Sa) f{a p + M p 5a)ti(M p )\0) - exp[ia p A (M p - I)6a/2h] (50) 

x(Sa\n(M p )\5a) . 
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Consider the expansion of the phase in ( |4q ) about the phase value R p characterising 
the period orbit. From (|4^), the difference S^/t between the phases characterising two 
nearby trajectories satisfies the following differential equation: 

S^t = — \5a t A a t — |a t A 5a t — \5a t A 5a t — 5H = —^^( a t A M t 5a) . (51) 

This equation must be solved subject to the boundary condition that J70 = 0. The 
solution is 

Sjt — — \ a t A MtSat + TjOtQ A 5a (52) 

For a periodic orbit, a t = olq = a p , so that in the neighbourhood of such an orbit 
(with M p = M t ) 

67 = -~a v A (M p - I)5a . (53) 
Expanding the phase in ji^ ) in the vicinity of a pi one now obtains 

R t - Pt q t /2+ PQ q /2 = R p - |[a p A (M p - T)6a] + 0(<5a 3 ) . (54) 

Taking everything together (note that all terms linear in 5a in the phase cancel) we 
obtain for the contribution Tr[L^] p of a periodic point a p to the trace of Ut- 

Tr[J7 t ] p ~ exp(ii? p /?i) Tr [K(M p )} . (55) 

It thus turns out that the stability property of the periodic orbit is encoded in the 
trace ofR.(M p ). Using ( |32"|) we have 

^.ap^/li + M) (56) 
v /|det(M p -/)| 

This contribution is summed over all periodic points and one obtains 

Tr £ ~ CXp(i ^ /n + i7r ^ /2 \ (57) 

periodic points 
of order n 

where p is now an integer index labeling the periodic points. The sum is over all 
periodic points, including primitive points and their repetitions. In the latter case, 
Tip — n/r p where r p is the order of repetition. With this result one finds for the 
density ( |l3|) of eigenphases 

d ,QS Ji Li >p n P exp[i(fl p - n p r p 8)/h + mv p /2] 

2 * n ^\det(M p -T)\ 

This is Gutzwiller's trace formula for chaotic maps, discussed in reference pj[] . The 
first term is often called Weyl contribution and describes the average density of states. 

4-2 Gutzwiller's trace formula for chaotic Hamiltonian flows 

In this section we consider chaotic Hamiltonian flows with d = 2 degrees of freedom 
and a time-independent Hamilton function H(q,p). As in the previous section, a 
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semiclassical approximation for the trace of the quantum-mechanical time evolution 
operator Ut will be derived, in order to find a semiclassical expression for the density 
of states. The trace Tr [Ut] will again be evaluated using coherent states, using (||). 
There are two contributions: a smooth part (Weyl term) due to very short orbits 
(corresponding to the first term in (|58|)) and oscillatory contributions from periodic 
orbits. 

We consider the smooth part first. For very short times, 

{a\U t \a) ~ (a\t(Sa t )\a) e"^"^ (59) 

where Sat — at — a. If the coherent states are those of the Hamiltonian H{a) = 
|a.a, then we may write 

(a\f(5a t )\a) = f(\6a t \/h) (60) 

where \a[ is the Euclidean norm of a, (and the function f(x) is a Gaussian, satisfying 
/(0) = 1). Thus the smooth part of the density is 

da f°° dt j[E-H(a)]t/h. 



(61) 

where f(y) is the Fourier transform of f{x). Now, noting that 

lim -f(x/e) = f(0)S(x) , (62) 

the contribution to the density of states from short times is 

(P(£)>~ J^sS[E-H( a )], (63) 

which is the well known Weyl estimate || . 

We now turn to the periodic-orbit contributions. The derivation of the previous 
section needs to be slightly modified since (as we shall shortly show) the stability 
matrix for a flow always has an eigenvalue unity (along the direction of the flow). 



Thus the representation (25) is singular. The solution to this problem is to treat 
longitudinal and transverse coordinates separately, as described in the following. 

It will be helpful to introduce a local canonical coordinate system in the vicinity 
of an unstable periodic orbit. The reference orbit is denoted by a p {r, rf), with r 
measuring the time taken a reach the point a p (r, rf) from some arbitrary starting 
point on the orbit at energy rj. The period of the orbit is T p (rf). Any point close to 
the orbit may be specified by giving its deviation from some point on the reference 
orbit: a — ct p (r, rf) + 5a, where 5a depends upon r and r/. For points close to the 
orbit at a reference energy 770, we will find it convenient to represent the deviations 
by a transformed set of local canonical coordinates, 5a' , such that 

5a = M'(r)5a' , 5a' = (St, SE, 5X, 5P) , (64) 
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where M' is a symplectic matrix depending upon r. This matrix is chosen so that 
the new coordinates 5a have the following interpretations. The energy 5E is de- 
fined by SE = H(a) — rjo. The time coordinate St measures the contribution to the 
displacement 5a from a vector parallel to the flow, namely JdH/da. The remain- 
ing coordinates 8(3 = (5X, SP) measure perpendicular displacements from the orbit; 
they are chosen such that Sf3 — corresponds to a point on the reference periodic 
orbit at energy t] + SE. The construction of the symplectic transformation M'(r) 
is described in an appendix. Note that this symplectic transformation is periodic: 
M'(T + T p ( m )) = M'(T). 

Given an initial deviation 8a from the reference trajectory a p (r, t]q) at time t = 0, 
there will be a corresponding deviation a t from a p (r + t, r) ) at a later time t. We 
write 

8a t = M'(t + T)8a' t , Sa' t = M' t '5a' , (65) 

where M' t ' is a monodromy matrix which describes the instability of the orbit in the 
local canonical coordinates da' . The construction of the Sa' coordinates implies that 
there is no coupling between the deviation in the (SE, St) subspace and the delta 5(3 
subspace. 

dt 

5a' = (St, Sr), 5(3) , 5a' t = (St + — Sr,, Srj, M ± S(3) , (66) 
so that the stability matrix M t " has the following form 

*-(*''£) <"> 

where M and M± are 2x2 symplectic matrices, with 

M\\ = ( I \), (68) 

and b = dt/dE, and M± is a 2(d— 1) x 2(d— 1) symplectic matrix mapping infinitesimal 
deviations in the transverse coordinates. In the fixed canonical coordinate system, 
the stability matrix describing the evolution of small deviations 5a from the reference 
trajectory is M t = M'(t + t)M£[M' '(t)}- 1 . 

Evaluation of the trace Tr[t/ t ] follows the approach used in section |4. l| , and re- 
quires the evaluation of {5a\TZ(M t )\5a) . The metaplectic operator lt(M t ) must be 
considered a a continuous function of t, defined in the extended domain. It may be 
written as a product, so that when t T p (r) ): 

Tr[H(M t )] = Tr[K(M' '(t + t))K(M?)[R.(M' '(t))}- 1 ] 

= (-l) s Tr[7t(Af'(r))^(Af t ")[^(M(r))]- 

- (-l) s Tr [£(<)] . (69) 

Here we have used the fact that M'(t) is per iodic, so that TZ(M'(t)) is periodic apart 



from a factor (— l) s discussed in section |3.3. 



The operator 1Z(M") may be written in the form 
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where TZ(M) and TZ(M±_) are operators acting upon Hilbert spaces corresponding to 
respectively one and d—1 degrees of freedom. These operators may be expressed as 
integrals of the form (p5|), constructed using sets of mutually commuting translation 
operators, respectively T(t, r/) and T({$). This implies that the matrix element of the 
metaplectic operator 1Z(M") is in the form of a product 

{6a\K(Ml')\6a) = {St, 6ij\K(M\\)\6t, Srf) (5 (3\R(M ±_)\5 (3) . (71) 

We now evaluate a contribution to Tr[£/ ( ] arising from a periodic orbit. The dominant 
contribution arises for period orbits of period T p = t, and we choose the energy r\ such 
that this condition is satisfied. Following the argument leading from ( fl8| ) to (p5|), the 
contribution to the trace is 

Tr[L> t ] p = J "dT J dS V J dS/3 (ot p + 5a\U t \a p + 5a) 



(2ttS) 2 

exp(LR p /7i + ins) 
2ttH 



dr / dSr) (5r) \H(M\ \)\8<q) Ti[n(M ± )} (72) 



where R p is the action of the periodic orbit, with the energy rj chosen such that 
Tp(ry) = t. Using Eqs. ([52]) and (|34|), and performing the integral over r, this 
evaluates to 

Tvrfn T p exp[iR p /h + mfa p + 2s)/2)] 

V 2mnb P v /|det(M±-/)|p 

where M± is evaluated for the periodic orbit with period T p = t, and /i p is an 
additional index chosen so that H(M±) is a continuous function of time along the 
periodic orbit. This expression is an explicit function of r\, but implicitly a function 
of t, since 77 satisfies t — T p (rj). 

Now the time integral over Tr[[/ t ] is performed, in order to evaluate the density 
of states. The contribution due to a single family of periodic orbits is 

f°° dt 

P(E) = J ^ eMiEt/h) Tr[U t } v . (74) 

This integral is approximated using the method of stationary phase. The time is 
written t = t a + St where to, the stationary phase point, satisfies 



dt 



(Et + R p ) = for t = t . (75) 



We will discuss the solution of this equation in detail. For each value of the time t 
over which we integrate, the integrand is evaluated at a different value of the energy, 
77, determined by the implicit relation t = T p (rj). The action S'p of the orbit and its 
period T p are also functions of the energy 77: 



S p (?]) = <j> p.dq 



T P = ^ • (76) 
(qr,p)=a p (r) or l 



The equation ([75]) for the stationary phase condition now becomes 

= ±{S p [ V (t)]-v(t)t + Et} = [T p (r ] )-i\^l + [E-r,(t)}. (77) 
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We have already required that T p (77) = t and we therefore satisfy the stationary phase 
condition by requiring, in addition, that E — i](t). When evaluating the integral (f74|), 
the value of the energy E therefore determines rj — E, which in turn determines the 
stationary phase point to — T p (rj). We note that this corresponds to a Legendrc 
transformation from R to S: 

S P = R P + Et, ^ = -E- (78) 

In order to perform the stationary phase integral, we require the second derivative of 
the action: we write 

Et + S p ( V (t)) - n t = S P {E) + ±A[i - T p (E)} 2 (79) 

where we find 



rl F 1 

*=-;=- = - 6-V (80) 



The Gaussian integration over 5t gives a factor of exp(i7r/2) Ak /27ri?i&p, and thus we 
obtain for the contribution of a periodic orbit to the density of states 

1 T p exp (i5 p /?i + i7rt/p/2) 



2nh '\det(M ± -I)\ 



where v p = fi p + 2s + 1 is the Maslov index of the periodic orbit. By summing over 
individual contributions from isolated, hyperbolic periodic orbits, labeled by an index 
p. Combining contributions from positive- and negative-time orbits one obtains: 

p{E) = (p(E)) + -^Re]T Tp expf^ + ^) . (82) 

nh P J\det(M ± -I)\ p ^ h 2 / 



This is Gutzwiller's trace formula Jl|, ||, |||. 

4-3 More general trace formulae for chaotic Hamiltonian flows 

The results of the previous section may be used to derive semiclassical approximations 
for more general traces, such as Tr[AC/t], where A is an observable with a classical 
limit A(a), compare Eq. (||). This trace may be written as 

Tr [AU t ] = J A(a) (a\U t \ a ) . (83) 

Using (^9|) one obtains for the smooth part of the density pa{E) [§o) 

( PA (E)) = J A(a) S[E - H(a)] . (84) 

As far as the periodic-orbit contributions are concerned, it is immediately obvious that 



all steps performed in section 4.2 are the same in this case, except for an additional 



factor of A(t, E) . One thus obtains for the weighted density |12| 

PA {E) = ( PA (E)) + ^-RcY: , A ' P expf ^ + ^) . (85) 
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where A p is the average of the observable A(a) over the p orbit: 

A p = f "dt A(a t ) . (86) 
Jo 

Similarly, a trace formula involving non-diagonal matrix elements may be derived, 
starting from the expression 

S A (t,s)=Tr[Au s AU t - s ]. (87) 

We may write 

5a( '' 5)= (2^/(5p^ )AK)(a| ^ |a> - (88) 
Using @, one has 

Jk eluJS J (2^p A{a) A{as) 5[E H{a)] ' (89) 

which is just the Fourier transform of the phase-space averaged correlation function, 
as derived in 0. Including the periodic-orbit contributions one obtains the following 
expression 

S A (E, M = (S A (E, M) + J_ Rc ^T S p (E,huj) ex Ji^ + ™A m 
where 

S p (E,hu;) = J — e lus J dt A(q t+S ,p t+S ) A(cc t ,p t ) . (91) 



5 Discussion 

We close with a couple of additional remarks. First, the above formulae were derived 
assuming that all periodic orbits are isolated and unstable (as is the case in hyperbolic 
systems). The above approach is very easily extended to deal with integrable systems 
where the dynamics is most conveniently expressed in angle and action variables 
a = (0,1). In these variables, the Hamilton function is a function of the action 
variables only, H{I). Thus 

I t = const. , e t = 9 + iv(I)t (92) 

with u>(I) = dH/dl. The for the typical case where the frequencies are not rationally 
related, the trajectories fill al dimensional tori in phase space. The periodic orbits only 
occur when the frequencies are rationally related, which occurs when the actions 
satisfy the implicit equation 

u[I(v)]t = 2icv . (93) 

where v — (v\, 1/2, ■ ■ ■) is a vector of integers. The trajectories labeled by the integer 
vectors v form d-parameter families (rational tori) which can be labeled by an initial 
angle 6q. It is thus clear that the leading-order contribution to the trace (||) comes 
from periodic orbits on rational tori. The derivation of the trace formula for integrable 
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systems [[TIJ proceeds exactly as before, using a 2d-dimensional generalization of the 
representation (|34|). 

Secondly, it is clear that the trace formulae Eqs. (0), (|S2"|), ( |85| ) and (|(j|) are not 
exact and at best asymptotically valid for small Ti. In order to obtain convergent 
expressions, the densities are often smoothed by convolution with a smoothing 

factor, i.e., 

Pe (E) = — / dE'p(E')cM-(E-E'f/2e 2 }. (94) 

This leads to a truncation of the periodic orbit sums in (|2|) , (|5|) and (|9(]) . The 
periodic-orbit formulae (^) and (^5|) are certainly valid provided the wave packet 
does not spread to the extent that the use of a linearised approximation to describe 
its deformation becomes invalid. If the system is hyperbolic, this criterion leads to the 
requirement that e > hX/loglSo/h), where A is an exponent describing the exponential 
separation of trajectories, and Sq a characteristic action scale for the system. In the 
case of (|9|), an additional smoothing with respect to the second argument of Sa is 
needed, of the same order of magnitude as e. 
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Appendix 

This appendix explains the construction of the coordinate system used in section 4.2 



It is constructed around a periodic orbit. The set of points on the periodic orbit is 
cx p (t, rj), where 77 is the energy and r is the time taken to reach the point starting 
from an arbitrary reference point. We assume that these reference points are chosen 
such that a(r, 77) is a smooth function of r. 

Local canonical coordinates (St, Srj, 5X, 8P) = (Sct'x, Sa' 2 , Sa' 3 , 8a' A ) = Set' are con- 
structed from the original coordinates (x,y,p x ,p y ) by obtaining a set of four vectors 
Vi, such that X)i v «^ a i i s the variation of (x,y,p x ,p y ) resulting from a variation of 
the new coordinates. The vector vi is the velocity vector of the Hamiltonian flow, 
vi = J(dH/dat). The vector V2 = da p (T,rj)/dr] satisfies VfJV2 = 1, so that the 
coordinates 5a[ = 5t and 6a' 2 — Sr) form a canonical pair. The remaining vectors 
V3 and V4 are constructed so that 7^ = vj Jvj satisfies 7^ = when i S {1,2} and 
j G {3,4}, and 734 = 1. These values of 7^ represent five linear constraints on the 
eight coefficients defining the vectors v 3 and V4. All of these vectors are functions of 
r. The matrix element M[Jr) of M'(t) is the i th component of the vector Vj. 
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Fig. 1: Illustrates the nearby orbit problem for a classical map. Shown is a periodic 
orbit ct p and a nearby orbit starting at a p + 5a and leading to a p + M p 5a in time t. 
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